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Let X = {X(t),t £ R*} be a random field with values in R d . For any finite Borel measure /i and 
analytic set E C R , the Hausdorff and packing dimensions of the image measure fj,x and image 
set X(E) are determined under certain mild conditions. These results are applicable to Gaussian 
random fields, self-similar stable random fields with stationary increments, real harmonizable 
fractional Levy fields and the Rosenblatt process. 
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1. Introduction 

Fractal dimensions such as Hausdorff dimension, box-counting dimension and packing 
dimension are very useful in characterizing roughness or irregularity of stochastic pro- 
cesses and random fields which, in turn, serve as stochastic models in various scientific 
areas including image processing, hydrology, geostatistics and spatial statistics. Many 
authors have studied the Hausdorff dimension and exact Hausdorff measure of the image 
sets of Markov processes and Gaussian random fields. We refer to Taylor (1986) and Xiao 
(2004) for extensive surveys on results and techniques for Markov processes, and to Adler 
(1981) and Kahane (1985) for results on Gaussian random fields. 

Let X = {X(t),t € 8 W } be a random field with values in R d , which will simply be 
called an (N, d) -random field. For any finite Borel measure \i on R^, the image measure 
of [i under X is defined by /ix Similarly, for every E C R , the image set is 

denoted by X(E) = {X(t),teE}c R d . This paper is concerned with the Hausdorff and 
packing dimensions of the image measures and image sets of random fields which are, in 
a certain sense, comparable to a self-similar process. Recall that X = {X(t),t G R^} is 
said to be H -self-similar if, for every constant c > 0, we have 

{X(ct),teM N } = {c H X(t),teR N } (l.i) 
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and X is said to have stationary increments if, for every /igl 



{X(t + h)- X(h),t e R N } = {X{t) - X(0),t G R N } 



(1.2) 



where = denotes equality of all finite-dimensional distributions. If X satisfies both (1.1) 
and (1.2), then it is called an H-SSSI random field. Samorodnitsky and Taqqu (1994) 
give a systematic account of self-similar stable processes. The main results of this paper 
show that the Hausdorff and packing dimensions of the images of an iZ-SSSI random 
field X are determined by the self-similarity index H and essentially do not depend on 
the distributions of X . 

An important example of an iZ-SSSI (N, d)-random field is fractional Brownian motion 
X = {X(t),t £ R N } of index H (0 < H < 1), which is a centered Gaussian random field 
with the covariance function E[Xi(t)Xj(a)] = ^6 id {\\s\\ 2H + \\t\\ 2H - \\t - s\\ 2H ) for all 
s,t £ 1^, where 8ij = 1 if i = j and Sij = otherwise. It is well known (see Kahane 
(1985), Chapter 18) that for every Borcl set E C R N , 



where dimn denotes the Hausdorff dimension. On the other hand, Talagrand and Xiao 
(1996) proved that, when TV > Hd, the packing dimension analog of (1.3) fails in general. 
Xiao (1997) proved that 



where dimp denotes packing dimension and Dim s E is the packing dimension profile of 
E defined by Falconer and Howroyd (1997) (see Section 2 for its definition). Results 
(1.3) and (1.4) show that there are significant differences between Hausdorff dimension 
and packing dimension, and both dimensions are needed for characterizing the fractal 
structures of X{E). 

There have been various efforts to extend (1.3) to other non-Markovian processes or 
random fields, but with only partial success; see Kono (1986), Lin and Xiao (1994), 
Benassi, Cohen and Istas (2003) and Xiao (2007). In order to establish a Hausdorff di- 
mension result similar to (1.3) for a random field X, it is standard to determine upper 
and lower bounds for dimn X (E) separately. While the capacity argument (based on 
Frostman's theorem) is useful for determining lower bounds, the methods based on the 
classical covering argument for establishing an upper bound for dimnX(E) are quite 
restrictive and usually require strong conditions to be imposed on X . As such, the afore- 
mentioned authors have only considered random fields which either satisfy a uniform 
Holder condition of appropriate order on compact sets or have at least the first moment. 
In particular, the existing methods are not enough, even for determining dimn A'QO, 1]^), 
when X = {X(t),t £ R N } is a general stable random field. 

Given a random field X = {X(t),t £ M. N } and a Borel set E C M. N , it is usually more 
difficult to determine the packing dimension of the image set X (E). Recently, Khosh- 
nevisan and Xiao (2008a) and Khoshnevisan, Schilling and Xiao (2009) have solved the 




a.s., 



(1.3) 



dim P X(E) = —Dim Hd E 



a.s., 



(1.4) 
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above problem when X = {X(t),t > 0} is a Levy process in M. d . However, their method 
depends crucially on the strong Markov property of Levy processes and cannot be applied 
directly to random fields. 

This paper is motivated by the need to develop methods for determining the Hausdorff 
and packing dimensions of the image measure fix and image set X(E) under minimal 
conditions on the random field X. By applying measure-theoretic methods and the theory 
of packing dimension profiles, we are able to solve the problems for the Hausdorff and 
packing dimensions of the image measure fix under mild conditions (namely, (CI) and 
(C2) in Section 3). The main results are Theorem 3.8 and Theorem 3.12. When X 
satisfies certain uniform Holder conditions, Theorems 3.8 and 3.12 can be applied directly 
to compute the Hausdorff and packing dimensions of X(E). More generally, we also 
provide a method for determining the Hausdorff dimension of X(E) under conditions 
(CI) and (C2) (see Theorem 4.9). However, we have not been able to solve the problem 
of determining dimp X(E) in general. 

The rest of this paper is organized as follows. In Section 2, we recall the definitions and 
some basic properties of Hausdorff dimension, packing dimension and packing dimension 
profiles of sets and Borcl measures. In Section 3, we determine the Hausdorff and packing 
dimensions of the image measure fix under general conditions (CI) and (C2). In Section 
4, we study the Hausdorff and packing dimensions of the image set X(E), where E C M. N 
is an analytic set (i.e., E is a continuous image of the Baire space N N or, equivalently, E 
is a continuous image of a Borel set). Section 5 contains applications of the theorems in 
Sections 3 and 4 to SSSI stable random fields, real harmonizable fractional Levy fields 
and the Rosenblatt process. 

Throughout this paper, we will use (x,y) to denote the inner product and j| • || to 
denote the Euclidean norm in R", no matter what the value of n is. For any s,t€ K n 
such that Sj < tj (J = 1, . . . , n), [s, t] = YVj=i i s j^j] ^ s called a closed interval. We will use 
K to denote an unspecified positive constant which may differ from line to line. Specific 
constants in Section i will be denoted by Kii, ifj,2, 

2. Preliminaries 

In this section, we recall briefly the definitions and some basic properties of Hausdorff 
dimension, packing dimension and packing dimension profiles. More information on Haus- 
dorff and packing dimensions can be found in Falconer (1990) and Mattila (1995). 

2.1. Hausdorff dimension of sets and measures 

For any a > 0, the a-dimensional Hausdorff measure of E C is defined by 




(2.1) 



Hausdorff and packing dimensions of the images 



929 



where B(x,r) = {y € M. N : \y — x\ < r}. The Hausdorff dimension of E is defined as 
diniH-E = inf{a > 0: s a -m{E) = 0}. For a finite Borcl measure p on M. N , its Hausdorff 
dimension is defined by diniR- p, = inf {dimn E: p(E) > and E C M. N is a Borel set} and 
its upper Hausdorff dimension is defined by dimjj p = inf{dimn E: p(M. N \E) = and E C 
is a Borcl set}. Hu and Taylor (1994) proved that 

dim H /x = sup(/3 > 0: limsup ^ B ^' r ^ = for /i-a.a. x e R N \, (2.2) 



r 



r 



dirnj^ p = inf <j /3 > 0: limsup ^ K y ^ > for /i-a.a. (2.3) 



The Hausdorff dimensions of an analytic set E C M. N and finite Borel measures on 
E are related by the following identity (which can be verified by (2.2) and Frostman's 
lemma) : 

dimH-E , = sup{dim H yu: /j£M+(£)}, (2.4) 

where JA^{E) denotes the family of all finite Borcl measures with compact support in 
E. 



2.2. Packing dimension of sets and measures 

Packing dimension was introduced by Tricot (1982) as a dual concept to Hausdorff dimen- 
sion and has become a useful tool for analyzing fractal sets and sample paths of stochastic 
processes; see Taylor and Tricot (1985), Taylor (1986), Talagrand and Xiao (1996), Fal- 
coner and Howroyd (1997), Howroyd (2001), Xiao (1997, 2004, 2009), Khoshnevisan and 
Xiao (2008a, 2008b), Khoshnevisan, Schilling and Xiao (2009) and the references therein 
for more information. 

For any a > 0, the a-dimensional packing measure of E C is defined as 

S a -p{E)=milY j 4>-P{E n ): Ec\jE n V 

^ n n ' 

where s a -P is the set function on subsets of R N defined by 

s a -P(E) = limsup< ^(2ri) a : B(xi,Ti) are disjoint, x% £ E,r{ <e>. 

The packing dimension of E is defined by dimp E = inf {a > 0: s a -p(E) =0}. It is well 
known that < dim H E < dim P E < N for every set E C R N . 

The packing dimension of a finite Borel measure pi on M. N is defined by dimp p = 
inf {dimp E: p(E) > and E C M. N is a Borel set} and the upper packing dimension of 
p is defined by dimp pL = infjdimp E: p(R N \E) = and E C is a Borel set}. In anal- 
ogy to (2.4), Falconer and Howroyd (1997) proved, for every analytic set E C M. N , that 



dimp E = supjdimp p: p 6 M.^(E)}. 



(2.5) 
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2.3. Packing dimension profiles 

Next, we recall some aspects of the packing dimension profiles of Falconer and Howroyd 
(1997) and Howroyd (2001). For a finite Borel measure fj, on and for any s > 0, let 

F?(x,r)= [ V s f— W(2/) 



be the potential with respect to the kernel 4>s(x) = min{l, ||a;|| _s },Va; £ M. N . 

Falconer and Howroyd (1997) defined the packing dimension profile and the upper 
packing dimension profile of /z as 

f F^(x r) 1 

Dim s a = sup<^ (3 > 0: liminf s V ' = for u-a.a. x € R w > (2.6) 

r— >o r' 



F^tx r) 

Dim* fi = mil P > 0: liminf s v ' ; > for /Lt-a.a. x £ R N [, (2.7) 



and 

n^-n* /, — inf J R n- i J±i 

respectively. Further, they showed that < Dim s /it < Dim* /i < s and, if s > N, then 

Dim s /i = dinip (j,, Dim* /i = dinip [i. (2-8) 

Motivated by (2.5), Falconer and Howroyd (1997) defined the s-dimensional packing 
dimension profile of E C by 

Dim s E = sup{Dim s /a: ^ e M+(.E)}. (2.9) 

It follows that 

0<Dim s £<s and Dim s E = dim P E if s > N. (2.10) 

By the above definition, it can be verified (see Falconer and Howroyd (1997), page 286) 
that for every Borel set E C R N with dinijj E = dinip E, we have 

Dim s £ = min{s,dimp£'}. (2.11) 

The following lemma is a consequence of Proposition 18 in Falconer and Howroyd 
(1997). 

Lemma 2.1. Let /j, be a finite Borel measure on M. N and E C R N be bounded and non- 
empty. Let a : R+ — > [0, N] be any one of the functions Dim s fi, Dim* /i or Dim s E in s. 
Then a(s) is non- decreasing and continuous. 
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3. Hausdorff and packing dimensions of the image 
measures 

Let X = {X(t),t G M. N } be an (N, d) -random field defined on some probability space 
(Cl,J-,F). We assume throughout this paper that X is separable (i.e., there exists a 
countable and dense set T* C M. N and a zero probability event To such that for every 
open set F cR N and closed set G C R d , the two events {lo: X(t, co)eG for all t G FnT*} 
and {w: X(t,ui) € G for all t <G F} differ from each other only by a subset of To; in this 
case, T* is called a separant for X) and (t,u>) >-> X{t,oS) is B{M. N ) x ^-measurable, where 
£(IR W ) is the Borcl cr-algebra of R N . 

For any Borel measure fi on R , the image measure /ix of yU under 1 1— >X(t) is 

Mx(-B) := yu{i € R^: G 5} for all Borel sets B C K d . 

In this section, we derive upper and lower bounds for the Hausdorff and packing dimen- 
sions of the image measures of X, which rely, respectively, on the following conditions 
(CI) and (C2). Analogous problems for the image set X{E) will be considered in Section 
4. 

(CI) There exist positive and finite constants H\ and /3 such that 

P{ sup \\X(s)-X(t)\\ >h Hl u\< K 3A u~P (3.1) 

4 s-t || <h ' 

for all t € M. N , h £ (0,ho) and u > uq, where ho, uq and 1^3,1 are positive con- 
stants. 

(C2) There exists a positive constant H2 such that for all s,t <E and r > 0, 

F{\\X(s)-X(t)\\ < \\s - t\\ H2 r} < K\ 2 mm{l 7 r d } 7 (3.2) 
where K3 2 > is a finite constant. 

Remark 3.1. Since (CI) and (C2) play essential roles in this paper, we will now make 
some relevant remarks about them. 

• Condition (CI) is a type of local maximal inequality and is easier to verify when 
the random field X has a certain approximate self-similarity. For example, if X 
is iJi-self-similar, then condition (CI) is satisfied whenever the tail probability of 
supi| s _ t || <1 H-X'(s) — X(£)|| decays no slower than a polynomial rate; see Proposition 
3.2 below and Section 5. It can also be verified directly for Gaussian or more general 
infinitely random fields by using large deviations techniques without appealing to 
self-similarity. 

• There may be different pairs of (Hi,/3) for which (CI) is satisfied. We note that the 
formulae for Hausdorff and packing dimensions of the images do not depend on the 
constant /3 > 0, it is sup{ffi: (CI) holds for some (Hi, [3)} that determines the best 
upper bounds for the Hausdorff and packing dimensions of the image measures. 
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• For every point t <G R , the local Holder exponent of X at t is defined as 

« (t >= S 4>0: ||io E^l.„ 

Condition (CI) and the Borel-Cantelli lemma imply that ax(t) > Hi almost surely 
(see (3.9) below). However, (CI) does not even imply sample path continuity of X. 

• In Section 4, the following, slightly weaker, form of condition (C2) will be sufficient: 

(C2') There exist positive constants H^ and K 3 2 such that (3.2) holds for all 
s,t€R N satisfying \\t- s\\ < I and r > 0. 

• Condition (C2) (or (C2')) is satisfied if, for all s,t £ M. N (or those satisfying 
||s — t\\ < 1), the random vector (X(s) — X(t))/\\s — t\\ H2 has a density function 
which is uniformly bounded in s and t. As shown by Proposition 3.3 below, (C2) is 
significantly weaker than the latter. 

The following proposition gives a simple sufficient condition for an SSSI process X = 
{X(t),t£M.} to satisfy condition (CI). More precise information can be obtained if 
further distributional properties of X are known; see Section 5. 

Proposition 3.2. Let X = {X(t),t £1} be a separable, H -SSSI process with values in 
M. d . If there exist positive constants (3 > and 3 such that Hf3 > 1 and 

¥{\\X(1)\\ >u} < J FT 3 ,3w~' 3 V«>1, (3.3) 
then there exists a positive constant 4 such that for all u>l, 

P( sup \\X(t)\\ >u\ < K 3i u-P. (3.4) 
L te[o.i] J 

In particular, condition (CI) is satisfied with Hi = H and the same j3 as in (3.3). 

Proof. Without loss of generality, we can assume d = 1 . Since the self-similarity index 
H > 0, we have X(0) = a.s. Let T* = {t n ,n > 0} be a separant for X = {X(t),te 
[0, 1]}. We assume that = to < ti < t-z < ■ ■ ■ < t n < ■ ■ -. For any n > 2, consider the 
random variables Yk (l-<k<n) defined by = X(tk) — X(tk-i). For 1 < i < j < n, let 
= Ei=i^fe- By the stationarity of increments and self-similarity of X and (3.3), we 
derive that for any u > 1, 



k—i 



>u\=F{ \X(1)\ > {f _ u -p [ < K 3 , 3 u-^ t] U-i) HP - (3-5) 



Thus, condition (3.4) of Theorem 3.2 of Moricz, Scrfling and Stout (1982) is satisfied with 
g(i,j) = tj — U-i, a = Hf3 and <j>(t) = t@ . It is easy to see that the non-negative function 
g(i,j) satisfies their condition (1.2) (i.e., g(i,j) < g(i,j + 1) and g(i,j) + g(j + l,k) < 
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Qg(i,k) for l<i<j<k<ri) with Q = 1. It therefore follows from Theorem 3.2 of 
Moricz, Serfling and Stout (1982) that there exists a constant ^3,4 (independent of n) 
such that for all u > 1 , 



p( max \X(tj)\ > u\=F< max 
ll<j'<« J J l<j'<n 



fc=l 



En 



>w> < (3.6) 



Letting n— > 00 yields (3.4), which, in turn, implies that (CI) holds for H±=H. □ 

Next, we provide a necessary and sufficient condition for an (TV, cQ-random field X = 
{X(t),t£ R N } to satisfy condition (C2) (or (C2')). For any r > 0, let 

Proposition 3.3. Let X = {X(t),t € R w } 6e a random field with values in M. d . Condi- 
tion (C2) (or (0,1')) then holds if and only if there exists a positive constant -K3.5 such 
that for all r > and all s,t£ M. N (or for those satisfying \\t — s\\ < 1) , 

( j, r (x)E(c i< ' x ' X ^- x ^^ t -^ H2 ) dx < R\ 5 min{l, r d }. (3.7) 

Remark 3-4- Since <fi r (x) = 0(||a;|| _2 ) as ||x|| — > 00, condition (3.7) is significantly 
weaker than assuming that (X(t) — X(s))/\\t — s\\ H2 has a bounded density and can 
be applied conveniently to SSSI processes. We mention that (3.7) is also weaker than the 
intcgrability condition in Assumption 1 on page 269 of Benassi, Cohen and Istas (2003). 
It can be shown that Theorem 2.1 in Benassi, Cohen and Istas (2003) still holds under 
(3.7) and their Assumption 2. 

Proof of Proposition 3.3. Note that for every r > 0, the function <p r (x) is non-negative 
and is in L 1 (IR d ). The Fourier transform of 4> r is 

M^) = Il( 1 -^) + Vz € M d . 

In the above, a + := max(a,0) for all o£K. Since z £ B(0,r) implies that 1 — (2r) -1 |2ij| > 

i, we have ±B(a, r ){ z ) £ 2 d (/) r (z) for all z € R d . Here, and in the sequel, 1a denotes the 
indicator function (or random variable) of the set (or event) A. By Fubini's theorem, we 
obtain 



»{\\X(s)-X(t)\\ < \\s-t\\ H2 r} <2 d 



X(t)-X(s) 



\\t-s\\ H * 

2 d I MxMc i{x ' X{t) - X{s))/llt - 4H2 )dx. 
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Hence, (3.7) implies condition (C2). On the other hand, we have <f) r (z) < l B / 2v /5r) ( z ) 
for all z £ R d . Consequently, 



E 



X(t)-X(s) 



< ¥{\\X(s) - X(t)\\ < 2Vd\\s - t\\ m -r}. 



\\t-s\\ H ' 

Therefore, condition (C2) implies (3.7). This completes the proof. □ 
3.1. Hausdorff dimensions of the image measures 

First, we consider the upper bounds for the Hausdorff dimensions of the image measure 

Proposition 3.5. Let X = {X(t),t £ M. N } be a random field with values in R d . If con- 
dition (CI) is satisfied, then for every finite Borel measure /i on R , 



d, -j-j- diuiH M 
d, jj- dim^j /it 



(3.8) 



Proof. Let A > 1//9 be a constant. For any fixed s £ R N and the sequence h n = 2 " 
(n > 1), it follows from condition (CI) that for all integers n > max{log(l//io): kT~2 u c/ A }i 



P{ sup \\X(t)-X(s)\\ >2- ffl ™(log2") A } <Kn~P x . 

Since X^^Li n ~^ X < oo, the Borel-Cantelli lemma implies that almost surely 

sup \\X(t)-X(s)\\ < (log2) x 2- Hin n x V?i>n , (3.9) 

||i-*||<2-» 

where no = no(w,s) depends on w and s. By Fubini's theorem, we derive that, for any 
finite Borel measure /i on R^, almost surely (3.9) holds for u-a.a. s £ M. N . 

We now fix an uj £ il such that (3.9) is valid for u-a.a. s £ M. N and prove that both 
inequalities in (3.8) hold. In the sequel, uj will be suppressed. 

To prove the first inequality in (3.8), since dimn"x < d holds trivially, we only need 
to prove that dimn l^x < ^ dimn Without loss of generality, we assume dimn u x > 
and take any 7 £ (0,dimn Mx)- Then, by (2.2), we have 

limsupr" 7 / l{|| J/ _ x ||< r } d/ix(y) = for /ix-a.a. x £ W l . (3.10) 

r-»0 Jm<i ~~ 
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limsupr~ 7 / ^{\\x(t)-x(s)\\<r} d//(i) = for /i-a.a. s <G R N . (3.11) 
r-X) Jm N 

Let us fix s £ such that both (3.9) and (3.11) hold. For any e > 0, we choose 
ni > no such that n x < 2 en for all n > n\. By (3.9), we can write 

/ ^{\\X(t)-X(s)\\<r}dfJ,(t) > V] / l{||jc(t)-X( s )|t<r}d^(t) 

Jl N n=ni J2-"- 1 <||t-s||<2-" 

(3-12) 

> / l { || t _ s ||< r i/(ff 1 - e) }d//(t). 

J||*-sj|<2-' l i 

Hence, we have 

l/|| t _ s ||< r i/(ffi- e )}d//(t) < / l{||x(t)^X( S )||<r}d/i(t) 

/TTD TV 

(3.13) 

^{llt-sll^rVCHi-s)} d/i(i). 



For the last integral, we have 



||t- s ||>2-"i 



limr-T/ l { || t _ fl ||< r i /( H 1 -.) } d^) = (3.14) 

'||t-s||>2-™i 



because the indicator function takes the value when r > is sufficiently small. 
It follows from (3.11), (3.13) and (3.14) that with r = p Hl ~ £ , 



lim sup p 



-(Hi-eb 



/ t {\\t-s\\< P } dfi(t) = limsupr 7 / lj|| t _ s ||< r i/(Hi-«n d(j,(t) 

(3.15) 



< limsupr 7 / l{||x(t)-X( s )||<r} d/i(t) = 0. 



We have thus proven that (3.15) holds almost surely for /i-a.a. s £ R^. This implies 
that dimH p > — e)7 almost surely. Since e > and 7 < dimn /^x are arbitrary, (3.8) 
follows. 

To prove the second inequality in (3.8), it is sufficient to show that dim^/ix < 
■^j-dimjj/i a.s. Let L) £ Q be fixed as above. We take an arbitrary (3 > dim^/i. By (2.3), 



we have 



limsup/O PI l{|| t _ s ||< p j dfj,(t) > for /i-a.a. s e R N . (3.16) 
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By using (3.12), we derive that for x = X(s), 

l{\\ y - x \\<r}dfix(y)> / t m _ sl i< r i/(H 1 -e) } diJ,(t) 



(3.17) 

1 {\\t-s\\<r 1 /' H i-'1} d(l(t). 

||t-s||>2 _ "i 

It follows from (3.17), (3.14) and (3.16) that 

for all s £ l" that satisfy (3.16). This implies that dim^jUx < PI {Hi — s) a.s. Letting 
£ J, and /3idimH/i yields the second inequality in (3.8). This completes the proof of 
Proposition 3.5. □ 



Remark 3.6. Note that in (3.8), the exceptional null probability events depend on fx. 
For several purposes, it is more useful to have a single exceptional null probability event 
fio such that, for all uj f2 , both inequalities in (3.8) hold simultaneously for all finite 
Borel measures fi on M. N . By slightly modifying the proof of Proposition 3.5 (see (3.12)), 
one can show that this is indeed true if, for every e > and every compact interval /, the 
sample function X(t) satisfies almost surely a uniform Holder condition of order Hi — e 
on I. 



Next, we show that condition (C2) determines lower bounds for the Hausdorff dimen- 
sions of the image measures of the random field X . 

Proposition 3.7. Let X = {X(t), t £ K^} be an (N, d)-random field satisfying condition 
(C2). Then, for every finite Borel measure li on M. N , 

diniH fix > mini d, — diurnal and dim^ fix > min i d, — dim^/jl a.s. 
I H 2 J I H 2 ) 

(3.19) 

Proof. In order to prove the first inequality in (3.19), we fix any constants < 7 < 7' < 
min{(i, -g- dimn fi}- Since dimu fx > "f' H 2l it follows from (2.2) that 

limsup ^^y^ =0 for ^-a.a. s£R Ar . (3.20) 

Let s £ 1* be a fixed point such that (3.20) holds. By (C2), we derive 

Efi X {B{X{s),r))= f F{\\X{t)-X{s)\\<r)fx{dt) 

jRN (3.21) 
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<K 3 MB( S ,r 1/H2 ))+K 3 , 2 f ( r ) fi(dt). 

J\\t-s\\>rV H * VII*- s \\ 2 J 

Let k be the image measure of fi under the mapping \\t — s\\ from M. N to R + . Then, 
by using an integration-by-parts formula and (3.20), we have 



t-s\\>r 



fr\ a f°° r d 

T^rMfl(-.P))dP (3-22) 

Jr 1 ' H 2 P 

< Kr 1 ' 

for all r > small enough, where the last inequality follows from (3.20) and the fact 
that i < d. Combining (3.21) and (3.22), wc sec that E/zjf (B(X(s),r)) < Kr' 1 ' for r > 
small. This, and the Markov inequality, imply that for all n large enough, 

P{nx(B{X(s),2- n )) > 2- n ~<) < K2- n{ ~<'-~<\ 

It follows from the Borel-Cantelli lemma that a.s. nx(B(X(s), 2"")) < 2~™ 7 for all n 
large enough. It should be clear the above implies that for all < 7 < min{d, dimu fi>}, 

Um sup f*x(B(x,r)) = Q for a x e Rd 

r^O ?" 7 

almost surely. Thus, dimnMA" > 7 a-S., and (3.19) follows from the arbitrariness of 7. 

To prove the second inequality in (3.19), let < 7 < 7' < min{d, dim H ^}. By (2.3), 
there exists a Borel set A C such that ^t(A) > and limsup^o r~ 7 H2 fi(B(s, r)) = 
for all seA The proof above shows that a.s. limsup r _ i . r _7 ^x(-B(a;, r)) = for all x G 
X(A). Since ^^(^(^l)) > a.s., we derive dim^ fix > 7 a.s. and the proof is completed. □ 

Combining Propositions 3.5 and 3.7, we have the following theorem, whose proof is 
omitted. 

Theorem 3.8. Let X = {X(t), t £ 1™} be an (N , d) -random field and let H be a positive 
constant. If, for every e > 0, X satisfies condition (CI) with H\ = H — s, some j3 = /3(e) > 
and (C2) with H2 = H + e, then for every finite Borel measure fi on l", 

diniH fix = mini d, — dimn fi \ and dim^ fix = mini d 7 — dim^ fi 1 a.s. (3.23) 
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3.2. Packing dimensions of the image measures 

We now study the problem of determining the packing dimensions dimp fix and dimp fix ■ 
The following upper bounds for the image measures are proved by Schilling and Xiao 
(2009). 

Proposition 3.9. Let X = {X(t),t £ M. N } be a random field with values in M. d . If con- 
dition (CI) is satisfied, then for every finite Borel measure /i on M. N , 

dimp fix < -rr- DiniHjd /i an d dimp fix < 77- Dim^ d fi a.s. (3.24) 
Hi Hi 



Similarly to Remark 3.6, we have the following. 



Remark 3.10. If, for every e > and every compact interval / C WL N , X(t) satisfies 
almost surely a uniform Holder condition of order Hi — e on /, then almost surely both 
inequalities in (3.24) hold for all finite Borel measures fi on R . 

For the lower bounds of packing dimensions, we have the following proposition. 

Proposition 3.11. Let X = {X(t),t £ M. N } be an (N,d)-random field satisfying condi- 
tion (C2). Then, for every finite Borel measure fi on 



siV 



dimp/ix > — DirnH^/x and dimp fix > 77- Dim# d fi a.s. (3.25) 
H-2 H2 

Proof. We only prove the first inequality in (3.25); the proof of the second one is similar. 
We may, and will, assume that Dim^d/i > 0. For fixed s £ R , Fubini's theorem implies 
that 

EF£ X (X(s),r) = f Emm{l,r d \\X(t)-X{s)\\- d }dfi(t). (3.26) 

J Si" 

The integrand in (3.26) can be written as 
Emm{l,r d \\X(t)-X(s)\\- d } 

(3.27) 

= P{\\X(t) - X(s)\\ < r}+E{r d \\X(t) - X(s)\\- d ■ l { ,| X (t)-x(.)||>r}}- 
By condition (C2), we obtain that for all s,t £ R w and r > 0, 

It-s 



F{\\X(t) - X(s)\\ <r}< K 3 , 2 mm\ 1, lu \. (3.28) 



Denote the distribution of X(t) — X(s) by r s t (-). Let v be the image measure of T Sit (-) 
under the mapping T:z H> ||z|| from M. d to K + . The last term in (3.27) can then be 
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f r d f°° r d 

1 d (3-29) 
J!—W{\\X{t)-X{s)\\< P }dp, 

r r 



<d 



where the last inequality follows from an integration-by-parts formula. 

By (3.28) and (3.29), we derive that the last term in (3.27) can be bounded by a 
constant multiple of 

r r d r P d 1 , 

1 11 " J (3.30) 

if r> \\t-s\\ H *, 

\t-s\\ H *' 



r d 


min 


pd+i 






{- 




I < 





if r< lit-. 



H 2 



It follows from (3.27), (3.28), (3.29) and (3.30) that for any < e < 1 and s,t<= R , 

Emin{l,r d \\X(t) - X(s)\\- d } < K 3fi wm{l , ^ ' _ e) }■ (3.31) 

For any 7 £ (0, Dim# 2t ;/j,), by Lemma 2.1, there exists £ > such that 7 < 
Dim ff2 ( d _ e ) /1. It follows from (2.6) that 

liminfr" 7/H2 / min/l,- r „ = , , , ldu(t)=0 for u-a.a. seR". (3.32) 

By (3.26), (3.31), (3.32) and Fatou's lemma, we have that for /i-a.a. s <E M. N , 
E^liminf r -i/ H ^F^ x (X(s),r)J 

<K 36 liminfr-T /ff2 / min/l,- * '.. „ . . , 1 d/x(t) = 

" ^ r^O 7 RJV \ p- S ||H 2 (rf-e)J ^ , 



(3.33) 

a— e 



By using Fubini's theorem again, we see that almost surely 

liminf r^ /H2 F^ x {X{s),r) = for /i-a.a. set*. 

Hence, dimp /Ltx > a.s. Since 7 can be arbitrarily close to DhiiH 2 d P, we obtain 
(3.25). □ 

The following is a direct consequence of Propositions 3.9 and 3.11. 
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Theorem 3.12. Let X = {X(t),t £ R N } be an (TV ,d)- random field and let H be a pos- 
itive constant. If, for every e > 0, X satisfies condition (CI) with Hi = H — e, some 
j3 = /3(e) > and (C2) with H2 = H + e, then for every finite Borel measure fi on R^, 

dimp fix = 77 Dini_ffd fi and dimp fix = — Dim* Hd fi a.s. (3.34) 



4. Hausdorff and packing dimensions of the image 
sets 

We now consider the Hausdorff and packing dimensions of the image set X(E). We will 
see that general lower bounds for dimH^(-E) and dimp X (E) can be derived from the 
results in Section 3 by using a measure theoretic method. For random fields which satisfy 
uniform Holder conditions on compact intervals, the upper bounds for dimn^(£') and 
dimpX(S) can also be easily obtained. However, it is difficult to obtain upper bounds 
for diniHW(-E') and dimpX(i?) under condition (CI) alone. We have only been able to 
provide a partial result on determining the upper bound for 6im^X(E). The analogous 
problem for dimpX(i?) remains open. 

We will need the following lemmas. Lemma 4.1 is from Lubin (1974), which is more 
general than Theorem 1.20 in Mattila (1995). 

Lemma 4.1. Let E CM. N be an analytic set and let f :M. N — > R d be a Borel function. If v 
is a finite Borel measure on M d with support in f(E), then v = fif for some fi G M+(.E). 

Lemma 4.2. Let E C be an analytic set. Then, for all Borel measurable functions 
f : R N — > M. d , we have 

diniH f{E) = sup{dim H fi f : ft £ M+(E)}, (4.1) 
dimpf(E) =sup{dimp/i/:/xe M+(E)}. (4.2) 

Proof. Denote the right-hand side of (4.1) by 7^. By (2.4), we get dimn f{E) > 773. 
Next, for any v £ M+(/(75)), Lemma 4.1 implies that v = fif for some fi £ M.^(E). This 
and (2.4) together imply dim.H/(-E) < je- Hence, (4.1) is proved. The proof of (4.2) is 
similar and is therefore omitted. □ 



We first consider the lower bounds for the Hausdorff and packing dimensions of X(E). 



Proposition 4.3. Let X = {X(t),t £ ~R N } be an (N,d)-random field that satisfies con- 
dition ( G2' ). Then, for every analytic set E C M. N , 

dimsX(E) > minid, -J-dimH-Ej and dim P X(E) > — Dim H2(J E a.s. (4.3) 
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Proof. Since both dimn and dimp are tr-stable (see Falconer (1990)), we may, and will, 
assume that the diameter of E is at most 1. Hence, condition (C2') will be enough to 
prove (4.3). 

Let us prove the first inequality in (4.3). It follows from (2.4) that for any < 7 < 
dirriH-E, there exists a fi € M+(£7) such that diniH/i > 7- By Proposition 3.7 (which 
holds for any finite Borel measure whose support has diameter < 1), we have dimn Hx > 
min{<i, -jj- dimn /1} a.s. This and (4.1) together imply that diniH^(-E) > min{<i, -^-7} 
a.s. Since 7 < dimn E is arbitrary, the desired inequality follows. 

Next, we prove the second inequality in (4.3). Note that for any < 7 < -g- Dim^ 2 d E, 
by (2.9), there exists a Borel measure fj, € M+(i?) such that H2J < Dim^^/i. It fol- 
lows from (3.25) that dimp /ix > 7 a.s. Hence, by Lemma 4.2, we have dimp X(E) > 7 
a.s., which, in turn, implies that dimpX(_E) > jj- Dim^^ E a.s. The proof is therefore 
completed. □ 

The following proposition gives upper bounds for dirriH^(-E) and dmipX(E). 

Proposition 4.4. Let X = {X(t),t G R N } be an (N ,d) -random field. If for every e>0, 
X satisfies a uniform Holder condition of order Hi — e on all compact intervals of R N 
almost surely, then, for all analytic sets E C K^, 

dim H ^(£) < mini d,^- dim E\ and dim P X{E) < Dim ffld E a.s. (4.4) 
I H 1 J Hi 

Proof. Both inequalities in (4.4) follow from Remarks 3.6, 3.10 and Lemma 4.2. □ 

Combining Propositions 4.3 and 4.4 yields the following theorem. 

Theorem 4.5. Let X = {X(t),t eR N } bean (N, d) -random field and let H e (0,1] be a 
constant. If, for every e > 0, X satisfies a uniform Holder condition of order H — e on 
all compact intervals ofM. N and condition (G2 1 ) with H2 = H + e, then, for all analytic 
sets EcR N , 

diniH^(-E) = mini d, —dimn e\ and dimp X{E) = — Dim^ E a.s. (4.5) 
I H J H 

It is often desirable to compute dimp X (I?) in terms of dimp E. The following is the 
packing dimension analog of (1.3). Note that if N > Hd, then the conclusion of Corollary 
4.6 does not hold in general; see Talagrand and Xiao (1996). In this sense, it is the best 
possible result of this kind. 

Corollary 4.6. Let X = {X(t),t e R N } and EdR N be as in Theorem 4.5. If either 
N < Hd or E satisfies dimn E — dimp E, then dimp X(E) = min{d, -k dimp E} a.s. 

Proof. If N < Hd, then (2.10) implies that for every analytic set E C M. N , Dim^ E = 
dimp E. Hence, Theorem 4.5 yields dimp X{E) = dimp E desired. On the other 
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hand, if an analytic set E C satisfies dim^E = dimpi?, then (2.11) implies that 
Dimnd E = mm{Hd, dimp E}. Hence, again, the conclusion follows from Theorem 4.5. □ 

Since many random fields do not have continuous sample functions and, even if they 
do, it is known that dimH^(-E') is not determined by the exponent of uniform modulus 
of continuity (a typical example being linear fractional stable motion - see Example 
5.4 below), there have been various efforts to remove the uniform Holder condition. 
However, except for Markov processes or random fields with certain Markov structure, 
no satisfactory method has been developed. The main difficulty lies in deriving a sharp 
upper bound for dimnX(E). 

In the following, we derive an upper bound for drain X(E) under condition (CI). This 
method is partially motivated by an argument in Schilling (1998) for Feller processes 
generated by pseudo-differential operators and, as far as we know, is more general than 
the existing methods in the literature. 

Lemma 4.7. Let X = {X(t),t G M. N } be a random field with values in M. d . If condition 
(CI) holds for H 1 > and /3 > 0, then, for allteR N , h>0 and 7 > 0, 

E(D(t, hye-^'V) < K 4tl h H ^^\ (4.6) 

where D(t, h) = sup|| s _ t || < ^ ||X(s) — -X"(i)|| and if 44 is a constant independent of t and 
h. 



Proof. We write 

/•oo 

E(D(t,hye- D ^' h) )= / u<- 1 c- u { 1 ~u)f'{D{t 1 h)>u}du 



(4.7) 

<K u 7 -V u (7-u)nfin{l,(/i- ffl u)~ /3 }du, 
Jo 

where the inequality follows from (CI). It is elementary to verify that, up to a constant, 
the last integral is bounded by 

h H i ,7 

v?- l du + h H ^ / uT- /3 - 1 ( 7 -u)du<if 4 i/i ffl( ' , ' A ' 3) . (4.8) 
Jh H i 

This proves (4.6). □ 

Proposition 4.8. Let X = {X(t),t £ R N } be a random field with values in M. d . Suppose 
that the sample function of X is a.s. bounded on all compact subsets ofR N . If condition 
(CI) holds for Hi > and (3 > 0, then, for every analytic set E C R that satisfies 
dim H E < PHi , 

dim H X(E) < mini d, — dim H E 1 a.s. (4.9) 
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Proof. Without loss of generality, we assume that E C [0,1]^. For any constant 7 € 
(diniH E, (3Hi), there exists a sequence of balls {B{t k , h k ), k > 1} such that 

00 

£;climsupB(t fc ,/i fc ) and ^(2/i fc ) 7 < 00. (4-10) 

fc^oo fc=1 

For a constant M > 0, let f2jvf = {uj: sup tS [ ^ N ||X(i)|| < M}. Since the sample func- 
tion of X(t) is almost surely bounded on [0, 1] , we have limjif-foo IP(^m ) = 1- Note that 
X(E) C limsupfc^oo B(X(t k ),D(t k ,h k )) and, by Lemma 4.7, (4.10) and the fact that 
7 < (iHi , we have 



^E(D(t k ,h k y^tn M )<e 2M ^E(D(t k ,h k )^ H ^e- D ^ h ^) 
fc=i fe=i 

4.11 

OO 

<C 2M ^4.1^^<0°- 



fe=l 

It follows from (4.11) that X)fcLi D(t k , hk) 1 ^ 1 < 00 almost surely on fijvf. This implies 
that diiriH -XX-E) < 7/^1 almost surely on Hm- Letting M — > 00 first and then 7 i dimn £7 
along the rational numbers proves (4.9). □ 

Putting Proposition 4.3 and Proposition 4.8 together, we derive the following theorem. 

Theorem 4.9. Let X = {X(t),t GM. N } be a random field with values in M. d whose sample 
function is a.s. bounded on all compact subsets of M. N . If there is a constant H > 
such that for every e > 0, X satisfies conditions (CI) with Hi = H — e and ( C2' ) with 
H2 = H + e, then for every analytic set E C K w that satisfies dimn E < /3H, 

dim H X(£;) =roinjd, dim H E I a.s. (4.12) 



5. Applications 

The general results in Sections 3 and 4 can be applied to wide classes of Gaussian or 
non-Gaussian random fields. Since the applications to Gaussian random fields can be 
carried out by extending Xiao (2007, 2009), we will focus on non-Gaussian random fields 
in this section. 



5.1. Self-similar stable random fields 

If X = {X(t),t € R+} is a stable Levy process in K d , the Hausdorff dimensions of its image 
sets have been well studied; see Taylor (1986) and Xiao (2004) for historical accounts. The 
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packing dimension results similar to those in Sections 3 and 4 have also been obtained 
by Khoshnevisan, Schilling and Xiao (2009) for Levy processes. In this subsection, we 
will only consider non-Markov stable processes and stable random fields. 

Let Xq — {Xo(t),t € ~R N } be an a-stable random field in R with the representation 



where M is a symmetric a-stable (SaS) random measure on a measurable space (T, .7-") 
with control measure m and f(t, •) : F — > R (t £ M. N ) is a family of functions on F satis- 
fying 



For any integer n > 1 and t%,...,t n £ R , the characteristic function of the joint distri- 
bution of Xo(ti), . . . ,Xo(t n ) is given by 



where £j GR (1 <j < n) and || • \\ a . m is the L a (F, J 7 , m)-norm (or quasi-norm if a < 1). 

The class of a-stable random fields with representation (5.1) is broad. In particular, if 
a random field Xq = {Xo(t),t £ ~BL N } is a-stable with a^l or symmetric a-stable, and 
is separable in probability (i.e., there is a countable subset To C M. N such that for every 
t £ M. N , there exists a sequence {tk} C To such that Xo(tfc) — > -Xo(f) in probability), then 
Xo has a representation (5.1); see Theorems 13.2.1 and 13.2.2 in Samorodnitsky and 
Taqqu (1994). 

For a separable a-stable random field in R given by (5.1), Rosinski and Samorodnitsky 
(1993) investigated the asymptotic behavior of P{sup tg [ ,i] N |-^o(*)| > u } as u ~ * 00 ( scc 
also Samorodnitsky and Taqqu (1994)). The following lemma is a consequence of their 
result. 

Lemma 5.1. Let Xq = {Xo(t),t £ R w } be a separable a-stable random field in R given 
in the form (5.1). Assume that Xq has a.s. bounded sample paths on [0,1]^. There then 
exists a positive and finite constant -ft^i, depending on a, f and m only, such that for 
all u > , 




(5.1) 






sup 




(5.2) 



Remark 5. 2. In the above lemma, it is crucial to assume that Xq has bounded sample 
paths on [0,1]^ almost surely. Otherwise, (5.2) may not hold, as shown by the linear 
fractional stable motion Xq with < a < 1 (see Example 5.4 below). 
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We define an a-stable random field X = {X(t),t £ R w } with values in R d by 

X(t) = (X 1 (t),...,X d (t)), (5.3) 



where X\ Xd are independent copies of Xq. 

The following result gives the Hausdorff and packing dimensions of the image measures 
of self-similar stable random fields. 



Theorem 5.3. Let X = {X(t),t £ R w } be a separable a-stable field with values in M. d 
defined by (5.3), where Xq is given in the form (5.1 ). Suppose that Xq is H-SSSI and its 
sample path is a.s. bounded on all compact subsets o/R . Then, for every finite Borel 



measure fi on R , 

diniH fix = mini d, — dimu^l and dimp fix = — Drniijd fi a.s. (5-4) 
I H J H 

Moreover, for every analytic set E C R that satisfies dimu E < aH , we have 
dimjj X(E) = min^c?, — dinifj E^ a.s. 

Proof. It follows from the self-similarity and Lemma 5.1 that X satisfies condition (CI) 
with H i = H and /3 = a. On the other hand, condition (C2) with H2 = H is satisfied 
because X is _ff-sclf-similar and has stationary increments, and the a-stable variable X(l) 
has a bounded continuous density function. Therefore, both equalities in (5.4) follow from 
Theorems 3.8 and 3.12. Finally, the last conclusion follows from Theorem 4.9. □ 



Next, we consider two important types of SSSI stable processes. 



Example 5-4 (Linear fractional stable motion). Let < a < 2 and H £ (0, 1) be 

given constants. We define an a-stable process Xq = {Xo(t),t £ M+} with values in R by 

X (t) = J h H (t,s)M a (ds), (5.5) 

where M a is a symmetric a-stable random measure on R with Lebesgue measure as its 
control measure and where 

M*> s) = a{(t - sf + - 1/a {-sf + - 1/a } + b{(t- ,s)5- 1/a - (-s)5~ 1/Q }. 

In the above, a, b £ R are constants with \a\ + \b\ ^ 0, t + = m&x{t, 0} and t- = max{— t, 0}. 
The a-stable process Xq is then if-self-similar with stationary increments, which is called 
an (a, H) -linear fractional stable motion. If H = — , then the integral in (5.5) is under- 
stood as aM([0,t}) if t > and as bM([t, 0]) if t < 0. Hence, Xq is an a-stable Levy 
process. 
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Maejima (1983) proved that if aH < 1, then Xo is a.s. unbounded on any interval of 
positive length. On the other hand, if aH > 1 (i.e., 1 < a < 2 and l/a< H <1), then 
Kolmogorov's continuity theorem implies that X$ is a.s. continuous. In the latter case, 
Takashima (1989) further studied the local and uniform Holder continuity of Xq. His 
Theorems 3.1 and 3.4 showed that the local Holder exponent of Xq equals H. However, 
the exponent of the uniform Holder continuity cannot be bigger than H — — . 

Now, let X = {X(t),t £ R+} be the (a, £f)-linear fractional stable motion with values 
in R d defined by (5.3). It follows from Theorem 5.3 that if aH > 1, then for every finite 
Borel measure /i on R + , 

diniH fJix = mini d, — dinin M 1 and dimp /ix = — Dim^ \i a.s., 
I H J H 

and for every analytic set E C R+, dimn^(-E) = min{<i, -g-diniH-E} a.s. Note that the 
above dimension results do not depend on the uniform Holder exponent of X. 

There are several ways to define linear fractional a-stable random fields; sec Kokoszka 
and Taqqu (1994). For example, for H e (0, 1) and a <G (0,2), define 

Z H {t)= ( {\\t- s\\ H - N ' a -\\s\\ H - N ' a )M a {ds) VteM N , (5.6) 

where M a is an SaS random measure on M. N with the iV-dimensional Lebesgue measure 
as its control measure. This is the stable analog of the ./V-parameter fractional Brownian 
motion. However, it follows from Theorem 10.2.3 in Samorodnitsky and Taqqu (1994) 
that, whenever N > 2, the sample paths of Z H are a.s. unbounded on any interval in 
R N . Thus, the results of this paper do not apply to Z H when N > 2. In general, little is 
known about the sample path properties of Z H . 



Example 5.5 (Harmonizable fractional stable motion). Given < a < 2 and H € 
(0, 1), the harmonizable fractional stable field Z H = {Z H (t),t g R^} with values in R is 
defined by 

Z H (t)=Re / rw ^^/^ (dA), (5.7) 

where M a is a complex-valued, rotationally invariant a-stable random measure on R N 
with the TV-dimensional Lebesgue measure as its control measure. It is easy to verify that 
the a-stable random field Z H is i/-self-similar with stationary increments. 

It follows from Theorem 10.4.2 in Samorodnitsky and Taqqu (1994) (which covers the 
case < a < 1) and Theorem 3 of Nolan (1989) (which covers 1 < a < 2) that Z H has 
continuous sample paths almost surely. Moreover, it can be proven that Z H satisfies the 
following uniform Holder continuity: for any compact interval I = [a, b] C M. N and any 
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£>0, 



\Z H It) - Z H (s)\ , \ 

hm sup n M /„ , i / ; = a.s. (5.8 



•i-t\\<h 



When N = 1, (5.8) is due to Kono and Maejima (1991). In general, (5.8) follows from 
the results in Bierme and Lacaux (2009) or Xiao (2010). Note that the Holder continuity 
of Z H is different from that of the linear fractional stable motions. 

Applying Theorem 4.5 to the harmonizable fractional stable motion in M. d defined as 
in (5.3), still denoted by Z H , we derive that for every analytic set E C R^, 

dim H Z H (E) = minjd, dim H £ j and dim P Z H (E) = Dim Hd E a.s. (5.9) 

Remark 5. 6. The results in this section are applicable to other self-similar stable ran- 
dom fields, including the Telecom process (Levy and Taqqu (2000), Pipiras and Taqqu 
(2000)), self-similar fields of Levy-Chcntsov type (Samorodnitsky and Taqqu (1994), 
Shieh (1996)) and the stable sheet (Ehm (1981)). We leave the details to interested 
readers. 



5.2. Real harmonizable fractional Levy motion 

We show that the results in Sections 3 and 4 can be applied to the real harmonizable 
fractional Levy motion (RHFLM) introduced by Benassi, Cohen and Istas (2002). To 
recall their definition, let v be a Borel measure on C which satisfies § c \z\ p v{dz) < oo 
for all p > 2. We assume that v is rotationally invariant. Hence, if P is the map z = 
pc 10 i v (&,p) € [0,27t) x R + , then the image measure of v under P can be written as 
vp{A9,Ap) = d0v p {Ap), where dd is the uniform measure on [0,27t) and v p is a Borel 
measure on R + . 

Let N(d£, dz) be a Poisson random measure on M. d x C with mean measure n(d£, dz) = 
E(N(d£,dz)) = d£v(dz) and let AT(d£,dz) = iV(d^,dz) - n(d£,dz) be the compensated 
Poisson measure. Then, according to Definition 2.3 in Benassi, Cohen and Istas (2002), a 
real harmonizable fractional Levy motion (without the Gaussian part) X 1 ^ = {X^(t),t £ 
R N } with index H € (0, 1) is defined by 

X" (t) = J RN ^2Re(^ ^H+N/l z ) ^( d ^ dz ) for all t G H^. (5.10) 

As shown by Benassi, Cohen and Istas (2002), X$ has stationary increments, as well 
as moments of all orders; it behaves locally like fractional Brownian motion, but at the 
large scale, it behaves like harmonizable fractional stable motion Z H in (5.7). Because 
of these multiscale properties, RHFLM's form a class of flexible stochastic models. 
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The following equation on characteristic functions of X^ was given by Benassi, Cohen 
and Istas (2002): for all integers n > 2, alH 1 , . . . , t" £ R N and all it 1 , . . . , «" £ R, 

EapfiW'lf^') =expf / [e/-«^-l-/ n (e,2)]dfi/(d2!)Y (5.11) 

\ ~1 / WxC / 

where 

/ n (e,z) = i2Re^^ ||g|| g + w / 2 J' 
In particular, for any s,t £ and ugR, (5.11) gives that 

Eexp m °,V » l ; Ucxp -27t/ ^ ^r^rM dg , (5.12) 



K « v n*-sinia H+jv/2 

where, for every - 0( a; ) is defined by ip(x) = f^°(l — cos(xp))i i p (dp). Note that 

the function ip is non-negative and continuous. Moreover, up to a constant, it is the 
characteristic exponent of the infinitely divisible law in C with Levy measure v. For 
the proof of Theorem 5.7, we will make use of the following fact: there exists a positive 
constant K such that 

ip(x) > K~ 1 x 2 / p 2 v p (&p) for aU x£ [0,1]. (5.13) 
Jo 

This is verified by using the inequality 1 — cos a; > K" 1 ! 2 for all x £ [0, 1]. 

Theorem 5.7. Let X ff = {X H (t),t £ M. N } be a separable real harmonizable fractional 
Levy field in M d defined by (5.3), where X^ is defined as in (5.10). Assume that ip 
satisfies the following condition: there exists a constant S £ (0,1] such that 

^^->a s for all a>l and x£R. (5.14) 
ip(x) 

Then, for every analytic set E C M. N , 

dmiH X H (E) = mini d, — dimn-E 1 and dimp X H (E) = — Dimnd E a.s. 
I H J H 

(5.15) 

Proof. It follows from Proposition 3.3 in Benassi, Cohen and Istas (2002) that for every 
e > 0, X H satisfies almost surely a uniform Holder condition of order H — e on all compact 
sets of M. N . Hence, the upper bounds in (5.15) follow from Proposition 4.4. 

In order to prove the desired lower bounds in (5.15), by Proposition 4.3, it suffices to 
show that X H satisfies condition (C2') with H-2 = H. This is done by showing that there 
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exists a positive function g <G such that for all s,t <S M. N satisfying \\s — t\\ < 1, wc 

have 

^ e i(u,X(t)-X(s))/\\t-s\\ H ^ < for aU u e R d_ ( 5 16 ) 

This and the Fourier inversion formula together imply that the density functions of 
X(t) - X(s)/\\t- s\\ H arc uniformly bounded for all s,teM. N satisfying ||s - t\\ < 1. 

Since the coordinate processes A^, . . . , A^ are independent copies of X^ , it is suffi- 
cient to prove (5.16) for d = l. Note that, by (5.16), we can take g(u) = 1 for all |u| < 1. 
For any u such that \u\ > 1, condition (5.14) implies that 



^U- S |iw^ 2 J de - XH A^liit-inieii^ 



£ll>ill'- 



d£ 

(5.17) 



>K\u\ / </> - — v/9 dg 



\\\t-4 H U\\ H+N/2 



where 7 > 1 is a constant whose value will be chosen later. 

By a change of variable £ 1— >• 77H i — s|| _1 , we see that the last integral becomes 



\\t-s\\ N / 2 (l- C os((t-s)/\\t-s\\,r)))\ d V 



s 



IN 



„||>7 V |M|™/2 J\\t 

(5.18) 

>K r (i-co S ((t-s)/\\t- s \\, V )r 
- 4n> 7 iwi 2H+Ar 

where the inequality follows from (5.13), and we have used the fact that \\t — s\\ < 1 and 
taken 7 large. The last integral is a constant because the Lebesgue measure is rotationally 
invariant. Thus, we have proven that for |u| > 1, 

Ec^m X "^-X° {s) ^ <exp(-K 5 , 2 \u\ s ). (5-19) 

Therefore, when d=l, (5.16) holds for the function g defined as g(u) = 1 if \u\ < 1 and 
g(u) = e~ Ks ' 2 \ u \ if |u| > 1. This completes the proof of Theorem 5.7. □ 



We mention that Benassi, Cohen and Istas (2004) have introduced another interest- 
ing class of fractional Levy fields, namely, the moving average fractional Levy fields 
(MAFLF). Similarly to the contrast between linear fractional stable motion and harmo- 
nizable fractional stable motion, many properties of MAFLF's are different from those 
of RHFLM's. For example, the exponent of the uniform modulus of continuity of an 
MAFLF is strictly smaller than its local Holder exponent. Nevertheless, we believe that 
the arguments in this paper are applicable to MAFLF's. This and some related problems 
will be dealt with elsewhere. 
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5.3. The Rosenblatt process 

Given an integer m > 2 and a constant k G (1/2 — l/(2m), 1/2), the Hcrmite process 
ym,« = {Y m ' K (t),tG K+} of order m is denned by 

Y m ^(t) = K 5 , 3 J^ y*p L (s-u j ) K + - 1 dsjdB(u 1 )---dB(u m ), (5.20) 

where > is a normalizing constant depending on m and n only and the integral 
L TO is the m-tuple Wiener-Ito integral with respect to the standard Brownian motion 
excluding the diagonals {ui = Uj}, i ^ j. The integral (5.20) is also well defined if m = 1; 
the process is a fractional Brownian motion for which the problem considered in this 
paper has been solved. 

The Hermite process Y m ' K is £f-SSSI and H = l + m/j-j £ (0, 1). It is a non-Gaussian 
process and often appears in non-central limit theorems for processes defined as integrals 
or partial sums of nonlinear functionals of stationary Gaussian sequences with long-range 
dependence; see Taqqu (1975, 1979), Dobrushin and Major (1979) and Major (1981). 

It follows from Theorem 6.3 of Taqqu (1979) that the Hermite process Y m - K has the 
following equivalent representation: 



/■' p it(u I +-+u m ) _ i ™ 

Y m > K (t)=K 5A — — -T\\u 3 \ K - 1 Z G {du 1 )---Z G (du m ), (5.21) 

J R m l{Ul + ■■ ■ +U m ) 



where > is a normalizing constant and Zq is a centered complex Gaussian random 
measure on K with Lebesgue measure as its control measure. 

Mori and Oodaira (1986) studied the functional laws of the iterated logarithms for the 
Hermite process Y m ' K . Lemma 5.8 follows from Lemma 6.3 in Mori and Oodaira (1986). 



Lemma 5.8. There exist positive constants and K^^, depending on m only, such 

that P{max te[M |r m ' K (t)| >u}< cxp(-A\ 6 w 2 /™) for all u > K 5 . 5 . 

Using Lemma 5.8, one can derive easily a uniform modulus of continuity for Y m ' K . 

Lemma 5.9. There exists a finite constant K§j such that for all constants < a < b < 
oo, 

\Y m - K (t + s) - Y m ' K (t)\ 
limsup sup sup -77- — , <K S>7 a.s., (5.22) 

Ho a<t<b-ho<s<h h H {logl/h) m ' 2 

where H = 1 + ma — y . 

Proof. For every t > and h > 0, the self-similarity of Y m ' K and Lemma 5.8 together 
imply that 

P{\Y m < K (t + h) - Y m ' K (h)\ > h H u] < exp(-A" 5 , 6 u 2 / m ). (5.23) 
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Hence, Y m ' K = {Y m ' K (t),t > 0} satisfies the conditions of Lemmas 2.1 and 2.2 in Csaki 
and Csorgo (1992) with cr{K) — h H and f3 = 2/m. Consequently, (5.22) follows directly 
from Theorem 3.1 in Csaki and Csorgo (1992). □ 

The case m = 2 has recently received considerable attention. The process Y 2 - K is called 
the Rosenblatt process by Taqqu (1975) (or fractional Rosenblatt motion by Pipiras 
(2004)). Its self-similarity index is given by H = 2k. This non-Gaussian process in many 
ways resembles fractional Brownian motion. For example, since H > 1/2, fractional noise 
of Y 2 - K exhibits long-range dependence. Besides its connections to non-central limit the- 
orems, the Rosenblatt process also appears in limit theorems for some quadratic forms of 
random variables with long-range dependence. Albin (1998a, 1998b) has discussed distri- 
butional properties and the extreme value theory of Y 2,K . In particular, Albin (1998b), 
Section 16, obtained sharp asymptotics on the tail probability of max tg r 0i i] Y 2 ' K (t). Pipi- 
ras (2004) established a wavelet-type expansion for the Rosenblatt process. Tudor (2008) 
has recently developed a stochastic calculus for Y 2 ' K based on both pathwise type calculus 
and Malliavin calculus. 

We now consider the Rosenblatt process X 2,K with values in M. d by letting its com- 
ponent processes be independent copies of Y 2,K . The following result determines the 
Hausdorff and packing dimensions of the image sets of X 2,K . 

Corollary 5.10. Let X 2 ' K = {X 2 - K (t) 7 t e K + } be a Rosenblatt process in R d as defined 
above. Then, for every analytic set E C K+ , we have 

dim H X 2 > K (E)=mm\d,-^-dirn H E\ and dim P X 2 ' K (E) = ^- Dim 2K( j E a.s. 
(_ 2k J 2k 

(5.24) 

Proof. By Lemma 5.9, for any e > 0, X 2 - K satisfies a uniform Holder condition of order 
H — e (where H = 2k) on all compact intervals in R + . On the other hand, it is known that 
the random variable Y 2 ' K (1) has a bounded and continuous density (see Davydov (1990) 
or Albin (1998a)). Thus, X 2 ' K also satisfies condition (C2) with H 2 = 2k. Therefore, the 
two equalities in (5.24) follow from Theorem 4.5. □ 
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